The role of symmetry on interface states in magnetic tunnel 

junctions 



C Uiberacker and P. M. Levy 
3^ ■ Department of Physics, New York University, 

^ '. 4 Washington Place, New York, NY 10003 



in 



C/2 



^ 



(N 
> 

in 
in 

o 
i> 

o 



c^ 



o 



% 



(February 1, 2008) 



Abstract 



When an electron tunnels from a metal into the barrier in a magnetic tunnel 
junction it has to cross the interface. Deep in the metal the eigenstates for 
the electron can be labelled by the point symmetry group of the bulk but 
C^ , around the interface this symmetry is reduced and one has to use linear com- 

binations of the bulk states to form the eigenstates labelled by the irreducible 
representations of the point symmetry group of the interface. In this way 
$Jh ' there can be states localized at the interface which control tunneling. The 

j^ . conclusions as to which are the dominant tunneling states are different from 

that conventionally found. 



Magnetic tunnel junctions (MTJ's) are layered structures of the form ferromagnetic 
metal/insulator/ferromagnetic metal, and provide a very interesting area of basic research. 
The tunneling currents of MTJ's are very much influenced by the electronic structure around 
the interface between the electrodes and the insulating barrier. For example by inserting a 
few layers of Cu (dusting layers) between one of the ferromagnetic Co electrodes and the 
AI2O3 barrier it was found recently |[l|J^ that the TMR ratio is suppressed exponentially. 
Furthermore, by putting IML Cr between Co and the same barrier the TMR effect nearly 
vanishes but when inserting 3-5 ML of Co between the Cr layer and the barrier the TMR was 
restored again |^. There have been theoretical attempts to model MTJ's with such dusting 
layers by using simple semiclassical models [^] and tight-binding schemes that gave some 
insight on the role of Fermi surface mismatch between different metals and resulting quantum 
well states by defining a tunneling density of states. When quantum mechanics is used to 
study electrons in crystals it follows from the symmetry of the Hamiltonian that the energies 
and eigenstates can be labelled by the irreducible representations (IR's) of the space group. 
Because of the strong screening of the Fermi sea of a metal electrons are not aware of the 
interface until they are within 3-4 monolayers of it; therefore electron states inside a metal 
far from the interface are labelled by the IR's of the point group in the bulk. Close to the 
interface the electron states should be classified according to the IR's of the point group of 
the interface which is a subgroup of the point group in the bulk. In addition, states that are 
localized at the surface may form, e.g., on transition-metals with (100) surfaces [||. These 
two features can be seen in Fig.(^ by comparing the density of states for minority electrons 
using as a basis the bulk eigenfunctions dz2,Pz,s (left figure) and the true eigenfunctions 
(right figure) at an Fe( 100) /vacuum interface. Three layers away from the interface the true 
eigenfunctions are nearly identical to the bulk eigenfunctions, but in a range of 2 layers 
around the interface dz2,Pz, s cannot be used to describe the electron; rather one has to use 
appropriate linear combinations of these states to form the true eigenstates. We also note in 
these figures that the DOS peaks in the region of the interface and decay exponentially into 
the barrier but also decay inside the metal indicating that the minority states are localized 
at the interface. In this brief report we outline the procedure for identifying which states 
are primarily responsible for tunneling in MTJ's. 

In order to analyze the electronic states in a MTJ one calculates the density of states 
(DOS) for it. Within a Green's function formalism and a single particle picture the DOS 
operator can be written as 

p{E; 6) = ^{G^- G^) = ^ (G{E - z6) - G{E + tS)) (1) 

2m zm ^ ' 

by making use of the operators G^^"^^ corresponding to the retarded and advanced Green's 
function. Here, i denotes the complex unit and the resolvent G{z) is defined by the equation 

{z - H)G{z) = i (2) 

with z being a complex number, H being the Hamiltonian and I being the unit operator. 
Using the eigenfunctions 1$^) of the Hamiltonian one can write the resolvent in the spectral 
representation and the DOS operator becomes 



Using the properties of |$j)($j| it follows that p{E;S) is normal and therefore has real 
positive definite diagonal elements in any basis. Especially, in a real space representation 
expressed in an angular momentum basis the matrix elements 

Pir;E;6) = g (^ E '^.,z.(r)$,V(r)^^— i^^^^ j Y,{r)Yl,{r) (4) 

are real and positive for each value of r. 

The MTJ we choose to investigate is Fe(100)/Vac/Fe(100) with an underlying b.c.c. lat- 
tice and we concentrated on the DOS around the surface in Fe(100)/Vac. For calculating the 
DOS we used the ab-initio spin-polarized scalar-relativistic Korringa-Kohn-Rostoker (KKR) 
multiple scattering formalism for layered systems |jll| together with the local spin density 
approximation (LSDA) and the Gunnarsson-Lundqvist exchange-correlation potential |]12| . 



To obtain short ranged structure constants, a screening formalism with a screening potential 
of 2.0 Ry was applied 0,113 and self-consistent potentials were evaluated within the atomic 
sphere approximation (ASA). First we calculated the potentials and exchange fields, where 
we used 16 points along a semi-circle in the complex energy plane with the points distributed 
on a Gaussian mesh and 45 k\\ points in the irreducible surface Brillouin zone (ISBZ). The 
vacuum was simulated by six layers of empty atomic spheres which have an underlying b.c.c. 
lattice structure. For calculating the DOS we choose an imaginary part of 5 = 5 x 10"'^ Ry 
and took 1830 ^n points in the ISBZ for the minority and majority spin-channel. 

In KKR for layered systems the arrangement of atoms within each layer is decomposed 
into Wigner-Seitz cells. An arbitrary point r^ is then uniquely described by the origin of cell 
n in layer p, R^, and the position r inside this cell and therefore r^ = R^ -|- r. In the ASA 
each cell is then replaced by a sphere of equal volume and the potential inside the cell is 
assumed to be spherically symmetric. It follows that the eigenf unctions of the Hamiltonian 
inside a given cell can be written as a series in angular momentum with expansion coefficients 
depending only on the distance from the origin of the cell. Then the Green's function can 
be written as 

LL' 

+ jrir)Zrir')eir-r')]}YLir)Yl,ir') , (5) 

by using L = [l, m) as a shorthand notation, in terms of the scattering path operator 
T^i'r'^'^ and the regular solutions Zl^'^{r) and the irregular solution Jl^'^{r) of the Schrodinger 
equation which all depend implicitly on E + i6 p!l|JT3[] . Using the fact that the Hamiltonian 
is hermitean it immediately follows that G{E — i6) = &{E + i6) and inserting this into 
Eqn.(|I|) yields 



p{rlT:i;E;6) = ^^{[G{r':i,rl-,E + z6)y-G{rlr':i-,E + z6)} 

= ^ E {[GT£''ir', r)r - G-rir, r')} yi(f)F£,(f') (6) 



LL' 

and the diagonal matrix elements are 



LU 
LU 

+ \m[Zr{r)Jr{r)]]YL{i)Yv{i) (7) 

It is very common to calculate the L-resolved integrated DOS (simply called L-resolved 
'density of states') which can be obtained from the previous equation as follows 

IYp(E;5)|L;n;p = /drp^^P(r;E;5) (8) 

The matrix of expansion coefficients in Eq. (J^) for a fixed cell n in layer p, which can be 
defined in terms of its components by 

p-'^{r-E-6) = {prv{r;E-6)] , (9) 

will be invariant under all operations of the point group. Due to the fact that the point 
group is a subgroup of the rotation-refiection group C(3) the basis functions of the IR's 
of the point group will be linear combinations of spherical harmonics. Moreover it can 
occur that in the truncated (l,m) basis used for calculations, which was up to / = 2 in 
our case, several spherical harmonics belong to the same IR and therefore p^''P{r; E; S) will 
have off-diagonal elements connecting these functions. By grouping these basis functions 
that belong to the same IR p"^'^(r; E; S) becomes block-diagonal with one block per IR. Now 
each block that contains more than one basis function has to be diagonalized to find the 
contributions p^'^{r] E; 6) belonging to the given IR, where j labels the eigenfunctions of H. 
If the quantization axis is chosen normal to the interface (z-direction) then at least all states 
{l,m = 0), which are functions of z only, will belong to the identity representation (they 
are invariant under all operations of the group) and the corresponding eigenfunctions will 
be linear combinations of these states. This is valid for all 2D points groups and is therefore 
very general. Due to the fact that p"'^(r; E; S) depends on r and p the transformation that 
diagonalize it will also depend on r and p and therefore the basis formed by the eigenfunctions 
('eigenbasis') will gradually change in a MTJ when going from the metal into the barrier. 

In Fe(001)/Vac the underlying crystal structure is b.c.c. and for the (001) surface the 2D 
lattice in the layers parallel to the surface is quadratic. The point group of the 2D lattice 
is therefore the group C^^ (Schonfiies notation) or 4mm (International notation) and it has 
5 irreducible representations Aj, ^ = 1, ...,5 where the first four are ID and the 5th is 2D 
PJ. For an angular momentum basis up to / = 2, A2 is not realized. Ai is the identity 
representation which is realized by s, Pz and ^^2 independently and we will concentrate on 
the mixing of these functions. 

In Figs.(0)-(H) we plot the L-resolved integrated DOS defined in Eq.(||) for the surface 
layer and the first and second vacuum layer against the energy for the functions belonging to 
Ai. Furthermore, we plot the total Ai contribution in the bulk in order to see if some states 
that occur at the surface do not occur in the bulk and therefore may be localized. Majority 
contributions are plotted with the correct sign and minority contribution with a reversed 
sign. The left graph shows the DOS in the (/, m) basis and the right graph shows the DOS 
in the eigenbasis and we will concentrate on the (/, m) basis first. Looking at the left part of 



Fig.(0) one can indeed see that at 1.6 eV below the Fermi Energy Ef in the majority DOS 
and around Ef in the minority DOS there are states that have no counterpart in the bulk. 
To find out if these states are localized we have a look at Fig.(|^) and Fig.(^ which show the 
DOS in the vacuum. Indeed, one finds a state at £"/-!. 6eV in the majority DOS and around 
Ef in the minority DOS which has a similar shape as the one at the surface layer, however 
we note it is largely of s-character whereas the state in Fig.@ was mainly of dz^ character. 
From that one could conclude that the state found at the surface of the metal does not 
continue into the vacuum and is therefore not localized. Moreover, in Fig.(^) it seems these 
contributions consist of two states with s and Pz character. However, these conclusions are 
wrong because when we look at the DOS in the eigenbasis (right part of Fig.(0)) we see 
that the peaks mentioned above result from one eigenstate that extends from the surface 
into the vacuum. If this were not the case there would be a state at the surface of Fe that 
would be damped almost completely at the first vacuum layer and one in the vacuum that 
is damped similarly at the surface of Fe, and at some point in between Fe and vacuum the 
DOS from these states must cross. This, however, is not possible because for all 2D point 
groups the eigenfunctions belonging to Ai span ID invariant spaces and therefore degenerate 
eigenvalues cannot occur in the Ai block of p; otherwise an invariant space of dimension 
> 2 would exist. To illustrate this fact we plotted in Fig.(|^) the component resolved DOS 
for the minority spin channel versus r by taking the diagonal elements of p"''P{r] E; 6) in the 
(/,m) basis (left graph) and in the eigenbasis (right graph), r = denotes the center of a 
cell in the surface layer and the vertical broken lines denotes the boundaries between layers 
when going from the surface into the vacuum. One can see that when transforming from 
the {l,m) basis to the eigenbasis any crossings of states occurring between the surface and 
the first vacuum layer are removed. The resulting eigenstates exist at the surface and in the 
vacuum. While we have focused on the minority states, the majority DOS versus r behaves 
qualitatively in the same way. 

The discontinuity of the eigenstates across the boundaries of layers in Fig.(^) comes from 
the fact that in the ASA the DOS is symmetric around the center of the cell in a layer and 
therefore is the same on the left and right boundary in Fig.(|^). By using a full potential 
calculation these discontinuities are removed, but otherwise it should quahtatively yield the 
same results. 

A typical wavelength of an electron at Ef is in the order of the layer spacing and it is 
only meaningful to present the average of the DOS over a given layer, therefore we plot 
the integrated DOS for each eigenstate on each layer in Fig.(|I|). Besides the facts discussed 
in the introduction one can see that different states within Ai have equal decay rates in 
vacuum. This is also observed for states belonging to A5 and it follows that each IR has 
one decay rate which then does not depend on the basis used. The different decay rates for 
the different IR's have already been pointed out and investigated by others |9|JT0[|. However, 
when they did their analysis they only took the trace over each block in J^^n drp"'^(r; E; S), 
obtaining the total contribution from each IR, and did not analyze the functions within a 
block to find the true eigenfunctions of the Hamiltonian (which are linear combinations of 
IR basis functions) that change when going from the electrode into the barrier. 

For completeness we want to mention that the states dx^-y-2 and dxy which belong to 
A3 and A4 will already be the eigenstates of the Hamiltonian but the 2D IR A5 is realized 
by two sets of basis functions, namely {d^z, dyz) and {px,Py)- Functions that belong to the 



same row of the same IR will mix and this applies to the pairs dxz,Px and dyz,Py Both 
pairs form the same pair of eigenstates and the two corresponding eigenvalues are twofold 
degenerate. The contribution from A3, . . . , A5 to the surface state found above is about an 
order of magnitude smaller than from the eigenf unctions belonging to Ai near the surface 
and becomes negligible in large barriers due to their larger decay. 

In summary we defined a density of states operator and its real space representation 
and showed that it is very useful for analyzing the electronic states in a magnetic tunnel 
junction. We showed that especially the point group of the layers parallel to the interface 
between the electrodes and the barrier and its irreducible representations can give some 
insight by making the DOS matrix defined in Eq.(|^) block-diagonal. But we stress that 
finding the eigenbasis within these blocks cannot be provided by symmetry but has to be 
found by solving a secular equation. Concentrating on the identity representation gives us a 
better understanding how localized states at (001) and (111) surfaces/interfaces in the cubic 
Bravais lattices can occur and that they will be of {s,Pz) — d^^ character independent of the 
actual material used for the electrode and the barrier (the z direction is chosen normal to 
the surface/interface). The general feature of s — d^i surface states in bcc (001) surfaces of 



different materials found by experiment has been mentioned in 0. Like others, see e.g. pO 
they use the bulk band structure and project it onto the face corresponding to the interface 
when discussing these states. This procedure implicitly uses a fixed {l,m) basis and not the 
eigenbasis, therefore it cannot explain how these states are formed. Finally we want to point 
out that states localized at the electrode/barrier interface play a crucial role in tunneling, 
and in the case of Fe/vacuum alter one's expectations about the spin polarization of the 
tunneling current. While one would argue on the basis of the itinerant (bulk) states in the 
junction that the current has a majority spin polarization, electrons can be scattered into 
the localized states at the interface and thereby reverse the spin-polarization of the tunneling 
current. This is reflected in the DOS at the interface around Ef, see Fig.(|^) . 

This work was supported by the Defense Advanced Research Projects Agency and the 
Office of Naval Research (Grant No. N00014-96- 1-1207 and Contract Nos. MDA972-96-C- 
0014, and MDA972-99-C-0009 ). 
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FIGURES 



Fe{100) minority-DOS at Ef, i: 1 ,SPD 



Fe(IOO) minority-DOS at Ef, i: 1 ,HYB 




FIG. 1. Ai contribution to the integrated DOS at the Fermi Energy plotted on each layer in 
the IR basis (left) and eigenbasis (right) for minority states. Layer denotes the surface layer and 
layers >0 denote vacuum. 
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FIG. 2. Ai contribution to the integrated density of states plotted versus energy for the surface 
layer of Fe(100)|Vac in the IR basis (left) and the eigenbasis (right). The minority DOS is plotted 
with reversed sign. 
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FIG. 3. Ai contribution to the integrated density of states plotted versus energy for the first 
vacuum layer of Fe(100)|Vac in the IR basis (left) and the eigenbasis (right). The minority DOS 
is plotted with reversed sign and 5 times enlarged. 
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FIG. 4. Ai contribution to the integrated density of states plotted versus energy for the second 
vacuum layer of Fe(100)|Vac in the IR basis (left) and the eigenbasis (right). The minority DOS 
is plotted with reversed sign and 5 times enlarged. 



Fe(IOO) minority-DOS at Ef, i= 1 ,SPD 




Fe(IOO) minority-DOS at Ef, i= 1 ,HYB 
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FIG. 5. Ai contribution to the minority DOS matrix p{r; E; 6) at the Fermi Energy plotted 
versus r in the IR basis (left) and the eigenbasis (right), r = denotes the center of the cell in the 
surface layer and r is measured in the direction to the vacuum. 
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